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ABSTRACT 

Aircraft  minimum  flying:  speed:,  as  determined: by  actual  flight  test, 
is  published  in   aircraft,  handbooks  for  pilot  guidance^  The  test  flight 
results  &re  used  ta  determine:  and  confirm  takeeeffand  landing  speeds, 
field  lengths,  Teft-hancb  portion  of  the  maneuvering  envelopes  (V-n 
diagram),  etc..  Qeterminatdan  of  the  absolute -minimum  flying  speed  of 
an  aircraft  an   the  other  hand:,.,  has  not  been  of- prime  importance  in 
flight  test- 

In.  the  present  analysis:  digital  simulation  allowed  the  systematic 
study  of  not  onTy  the  minimum  flying  speed. assdefined  by  Federal  Aviation 
Regulations  but  aTso  the  absxxlute  minimum  flying:speed  attainable  in 
steady,  unaccelerated  flights.  The  study  included  such  effects  as 
deceleration  rate,,  rate  a-f  change  of  elevator  angle,:  .aircraft  weight 
and  pitch  moment  of  inertiat.. 

It  was  found  far  an  assumed  light-weight: fighter  aircraft  that  the 
absolute  minimum  fTying  speed: was  approximately  20  knots  less  than  the 
FAR  minimum  fTying  speed..  Moreover  the  FAR: minimum  fTying  speeds 
tended  to  be  quite  sensitive  to  rate  of  change  of  elevator  angle. 
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symbols: 

2" 

a  normal  flight  path  acrBileration,.  ft/see 

nfp 

21 

B  pitch  moment  af  inertxa,.  srlug-ft: 

C  location  of  the  aircraft  center  of"  mass: 

c~  mean,  aerodynamic  chord. 

CD  drag  coefficient,.  drag/qS 

C.  lift  coefficient,.  Tift/qS 

C.  maximum  IF  ft  coefficient 
max 

C/  apparent  Tift  coefficient,.  weight/qSI 

C.'  maximum  apparent  Tift  coefficient. 
max 

C.  (o)  static  Tift  coefficient  as  a-  function  of  "angle  of -"attack 

C,  rate  of  change  of  Tift  coefficient. for  ae change: in- 

6  elevator  angle 

C.  .  rate  of  change  of  Tift  coefficient,  for  a: change  in 

e  pitch  velocity. 

Cm  pitching  moment  coefficient.,,  moment/ qScT 

Cw(a)  static  pitching  moment  coefficient  as:  a:  function  of 

1  angle  of  attack  ' 

8C.  ,. 

C„.  angle  of  attack  damping  derivative,  — (rad~   ) 

a  3[aC 
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Cj^.  pitching  moment  coefficient  derivative  for 

6„  elevator  deflection 
e 

CM  pitch  damping  derivative,  : —  (rad.  ) 

q  lee 

3\2y 

Cj  thrust  coefficient,  thrust/qS 

C„  longitudinal  aerodynamic  force,  coefficient, .body,  axis 

C  vertical  aerodynamic  force  coefficient,  body  axis: 


r 


D  jproduct  erf"  trrertrcE  jfjyzdm,.  slug-ft 

6  elevator  angles,,  deg;-. 

• 

6  rate  af  change  of  elevator  angle,  deg/sec: 

et  local"  truncat-iiair  error- for  Runge-Kutta:  algorithm 

F  product  of"  inertia  Jxy  dm,  slug-ft  ' 

F  resultant  external"  forces  vector  relative;to:aircraft 

cerrter  of  mass 

6  resultant  external"  moment,  vector  about: the: aircraft  center 

of  mass 

2' 
g  scalar  accel eration"  due:  to  gravity,  ft/sec 

F  angular  momentum  vector- relative  to:  aircraft  .center  of 

mass 

[L,M,N]  scaTar  components:  off (35;  rolling,  pitching  and  yawing 

moments,.  ft-Tb... 

[P,Q,R]  scaTar  components  of" W;:  rolling,  pitching:and  yawing 

veT.oci  ties ,.  rad/  sec . 

Q  rate  of  change  of"  piitching  velocity 

_  22 

q  dynamic  pressure,.  Itt/ft 

_  2 

q  dynamic  pressure  at  minimum  flying  speed/stall   speed,  lb/ft 

2Z 
S  wing  area,,  ft 

T  thrust,.  Tb. 

[u,v,w]    scalar  components  of  V  ,  relative  to  body  axis  system, 
ft/sec.  a 

2 
U        rate  of  change  of  V,.  ft/sec 

Y  rate  of  change  of  V  ,  as  determined  by  FAR  method 

Y  rate  of  change  of  V  ,  as  determined  by  altitude  break 
n        method  c 

Yp        magnitude  of  resultant  linear  velocity  of  aircraft  center 
mass*   ft/sec. 

Y  resultant  velocity  vector  of  aircraft :center  of  mass, 
g        ft/ sec 
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V  nrnrrfmum  fTyfmp  speed:  ass  determined: by.  CI  method;  .ft/sec. 
cmin 

V  ,  nriinrf'mum  flying:  speed:  ass  determined:  by  CI  method;  .ft/sec. 
nin 

V.  minimum  fTying.  speed:  ass  determined:  by.  altitude  .break 

nin  method,,  ft/ sec.. 

V  .  minimum  flying  speed:  as:  determined- by  constant  deceleration 
mn  irate*  ft/ sec 

W  aircraft  weight,.  The 

—  2 

W  irate  of  change  ofllv  scalar- componentsofv^,  .ft/sec   . 

[x,y,z]  cartesian  coordinate  axes:  notation  for  body  axis  :system, 

enrgin  located  at.  aircraft,  center  of- mass: 

[x'.y'jZ1']    cartesian  coordinate  axes- notation  for  stability  axes 
system,,  orgin  located:  at.  aircraft: center  of  mass  : 

[X,Y,Z]  components  of  resultant:  aerodynamic: force. acting, .on 

aircraft,,  lb.. 

y  +,  ,  variable  in  Runge-Kutta  algorithm  obtained  by  intergrating 

n     *  between  two  points:,.  X     and:  X     ,,  with  step  sizeh]. 

'  +1  7  variable  in  Runge-Kutta:  algorithm  obtained.by  intergrating 

n     *  between:  twer  points,.  X    and  X     ,-,-with  step  size. h9  where 

h2  =  hT/Z  n  nl  £" 

o  angle  of  attack,,  dec. 

o  rate  erf  change  of  angle,  of-  attack , . deg/sec . , 

e  afrcraft  pitch  angle,,  angle  between  horizontal   reference 

and  the  X  axis  of  the  body  axis  system,  deg. 

6  rate  of  change  of  pitch  angle,  deg/sec 

3 
p  dtir  density,  slugs/ ft 

[f,Q,$]  Euler  angles,   rad. 

to  anguTar  velocity  vector  of  the  aircraft,  rad/sec 
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E..  INTKOPUCTION' 

Minimum  flying  speed  and  or  stalling: speed: is: defined  by  Federal 
Aviation  Regulations  (FAR)  pfef"..  1*1  as:  the: minimum  steady  speed  :at" 
which  the  aircraft  remains  control  Table..  The:  NayaT  Test  Pilot  .School 
Flight  Test  Manual  [Ref..  Zl  defines:  minimum  fiying_ speed  as  : the  .minimum 
steady  airspeed  attainable  in  unacceTerated: flight: or  the .minimum 
usable  airspeed.  The  actual'  determi  nation  of'minimum  flying  .speed  is 
always  accomplished  through  flight  tests:. .  Methods  for  theoretical 
predictions  are  empirrcaT,.  based  upon  windrtunnel' results,  experience, 
etc. 

FAR  specified  that  stalls  would  be:  demonstrated: by.  trimming  :the 
aircraft  at  1:30  percent  of^  the  estimated  minimum  flying  speed  :and . 
decelerating  at  a  constant  rata- until"  a,  minimum  speed. was  .obtained. 
The  certified  stall  speed  would  correspond,  to  a: deceleration  rate  of 
1  knot/sec.  To  faciTftate  the  test  pilot' s:  task,. Ref ;  1  "and. 2 Idefined 
aircraft  characteristics  and  parameters.  indicative:of -"minimum  flying 
speed.  Typical  indications  were:  longitudinal,  directional  or  lateral 
divergence,  excessive  altitude  loss,  Toss:  of  control/effectiveness,  etc, 
Data  analysis  methods  and  restrictions  on  demonstration  technique, 
other  than  rate  of  change  of  airspeed,  were:  not: specified  by  Ref.  1. 

Additional  flight  test  techniques  have:  been  proposed  and  evaluated 
by  the  National  Aeronautics  and  Space  Administration  [Refs.3."and  4]. 
These  studies,  through  flight  tests,  evaluated  three  demonstration 
techniques:  FAR,.  T-g  break  and  constant-rate-of-elimb.  Three  methods 
of  data-analysis  were  used  with  the  FAR~ demonstration  technique  .to 
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determine  minimum  flying  speed::    constant:  dece-i) eraii on-, .  C?  and:C^  methods. 
The  minimum  flying  speeds  ahtained.  by  the:  three:  tect"im'ques:were-:corrected 
to  a  deceleration  rate  of  1  knot/ sec  ifr  accordance:  with  Ref; .  1". . 

The  study  described  herein  applied  numerical"  methods: to: simulate ■_ 
flight  test  data  in  an  attempt  to  define  an  absolute: mini mum' flying; 
speed.     Five  data-analysis  methods  were  used:..    The:  speed-was: an- absolute . 
minimum  obtainable  for  steady,,  un accelerated  flight". when- the; FAR: 
correction  to  a  1  knot/sec.  deceleration  rate:  was:  disallowed. .   The-: 
absolute  minimum  flying  speeds  were  compared  with-  FAR:  certified.minimum 
flying  speeds,  i.e.,  corrected  to  11  knot/ see-,  deceleration  rate i . 

"Flight  test  data"  for  the  F-34A,.  a  small"  s.ingl"e-engine£jet:ffghter, 
v/as  obtained  from  a  computer  program  of"  the:  non-linear" equations: of 
aircraft  motion.     FORTRAN  IV  language  and  the:  Naval'  Pbstgraduate-:School 
IBM  360  digital  computer  were  utilized..    The:  computer  program  provided: 
numerous  options   for  trimming  and  flying  the  aircraft,,  including.: 
variations  of  thrust,  weight,,  rate  of  change  of- elevator  angle -and: 
flight  orientation,  i.e..,.   climbing,,  level   or  descending-flight: 

For  this  study,  the  aircraft  was:  trimmed  irr  steady, .descending 
flight  at  130  percent  of  the  estimated  minimum  flying  speed  and  decelera- 
ted by  decreasing  the  elevator  angle  at  a  constant  rate.     Two  values . 
for  thrust  were  used  to  trim  and  stall   the  aircraft:     T/W  s.0.02.and 
T/W  =  0.11. 

The  effects  on  minimum  flying  speed  of  thrust,. aircraft  weight, 
rate  of  change  of  elevator  angle  and  aircraft  moment: of  inertia:about 
the  lateral   axis  were  also  studied. 
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The  purpose  af  thts  study,  was;  toe  invest! gate: the :•  relationships  of 
the  aircraft  minimum  flying  speeds:  as:  determined:  by.  the -.different  data- 
analysis  methods.. 
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IE..  EXPERIMENTAL  METHODS: 

A.  BASIC  APPROACH 

The  study  of  minimum  flying  speed  was:  restricted:  to:  aircraft  .motion 
in  the  XZ  plane:     i.e..*  to  aircraft:  motion  a:l~ong:  ttie:  X  and:Zlaxes:and 
pitching  moments  about  the  Y  axis..     I!t:  was:  assumed:  that: no:  adverse: 
lateral  or  directional  traits  occurred:  while:  decelerating: to: the: 
minimum  flying  speed. 

Non-linear  equations  of  aircraft  motion  with-  three:- degreesfef-:- 
freedom  were  derived  fron  EuTer's  equations:  of 'motion,.  Appendix;  A.  . 
Extensive  use  was  made  erf  Etkin  [Ref..  53I--    Thee  equations.  wereeset:up. 
for  input  of  non-linear  aerodynamic  data,,  by  means:  of  "a:  table:  look-up 
procedure. 

A  computer  program  was  written  trr  FORTRAN  IV.  language.:  fbr-the-iNctval  * 

Postgraduate  School   IBM  360  digital"  computer,.  Appendix  E£.   The^program 

employed  a  fourth-order  Runge-Kutta  algorithm  in solving  the:equations 

of  motion.     Real-time  was  used  in  evaluating  the: aircraft: motion  from 

approximately  1.3  V    •         to  V    .  Typical' time:  histories  of  the 

aircraft  parameters  were  plotted  in  Figure  T.     Figure. 2. contains  summary 

plots  of  V   .     as  a  function  of  deceleration  rate  for  differentiates 
min 

of  change  of  elevator  angle.     The  thrust  to  weight  ration   (T/W)  was  0.02 

B.  TEST  AIRCRAFT 

The  F-94A,  a  small   single-engine  jet.  fighter,,  was:  utilized  as  the . 
test  aircraft. 

Reference  aerodynamic  data  were  obtained  from  Blakelock   [Ref; "  .53 "in 
the  form  of  linear  stability  derivatives  for  all   terms  except: those 
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Figure  1.     Typical  Time  Histories  of  Aircraft  Parameters 
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depending  upon  angle  of  attack. .    The;  variation  of 'tbe;complete-:aircraft 
C.  ,  CD    and  C«»  was  estimated:  from  thee  avaiTabl  e  data: for  a^center  of 
gravity  located  at  the  Q..Z52  Z  (wing:  mean  aerodynamic: chord)  station 
and  an  airplane  in  the  landing  configuration.     Although  the^variations 
of  C,  ,  CD  and  CM  with  angle  of  attack,,  a,  were  estimatedfor  zero 

elevator  deflection*  the  use  of  the:  linear  stability  derivatives  such 

C       C 

as  L.  and  M~  allowed:  trimming;  the:  aircraft  to  any  initial  flight 

e       e 

conditions.    The  C|;   vs..  at  curve  for  6^  =0  was  selected:for  a:trim 

L  ee 

C.         of  about  IAS  with  the.  stalT  break  being  very  gradual,  .which  is 

max 
indicative  of  a  progressive  flow  separation  starting  at~the:wing  trailing 

edge. 

Aircraft  reference  data  and  physical"  characteristics:  are -listed  in 
Table  I. 

The  choice  of  aircraft  fTrght  traits:  are  arbitrary  and  :may  seem 
restrictive  to  a  person  who  is  solely  concerned  about;a: particular 
aircraft  by  virtue  of  personal  experiences.  However,  it.should  be 
recognized  that  the  methods  and  qualities  described,  in-  this  analysis 
are  typical  and  may  be  modified  to*  a  particular  application  by  changing 
the  input  data. 

C.  FLIGHT  PROCEDURES 

All  "test  flights"  were  conducted  with  the  F-94A  configurated  for  a 
landing  approach.  The  aircraft  was  trimmed  for  steady  descending  flight 
at  130  percent  of  the  estimated  minimum  flying  speed.  Aircraft  thrust, 
weight  and  moment  of  inertia  about  the  lateral  axis. were  the  only 
physical  characteristics  changed  during  the  "test  flights". 
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TABLE   I 

a-    PbysfcaV  Data 

Parameter  Value 

$  239:ftr 

c  6t4-ft": 

W,  12;359.'lb 

E  26i545islug-ft 

(L,   (Oinear  region)  5;2Z~rad~ 


a 


X 


a 


\ 


St. 

c 


2 


0:43  2rad'] 


(L.,  (O.ifiear  reg.ibir)  -0:406  rad" 


-QJSS.rad"1 


-4r268  rad"1 


-8.165  rad"1 


q 

C  25;2%  c 
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TABLE  E  ((Cont'd) 


b.     Tabular  Data  (complete  aircraft,  6'    =  03 

e 


(rad) 

^ 

c^ 

ck 

0.10472 

CL552 

-<L0D27" 

0:0940 

0.13963 

Q..738 

-cloito: 

o:ii5o 

0.17453 

CL92T. 

-0:0315: 

o:i4oo 

0.20944 

Lags 

-0. 0445: 

0:1695 

0.24435 

T:..Z32 

-010585: 

012030 

0.27925 

1L346: 

-<r.074o: 

0:2490 

0.31416 

I. .430 

-0:09m 

0:3120 

0.34907 

I -484 

-(LMW 

0'.3640 

0.38397 

T.5T:3 

-0U32D: 

014050 

0.41888 

T..52H 

-a.iB2D: 

0.4550 

0.45379 

T.sao 

-0U720" 

0:4930 

0.48869 

I..462 

-(L2D75: 

0C5350 

0.52360 

T.4T.5 

-a.  24oo: 

0.5720 

0.55851 

1.360 

-a.  2795: 

o:6ioo 
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Deceleration  from  the  trtnr  exmrtttiarr  was:  accompli  shed:  by.  decreasing 
the  elevator  angle,  to  simulate,  aft  cnntral  stick-  movement..  Available: 
time  histories,  [Ref.  41,.  indicated,  elevator-  angle:  to:  be: a: near  'linear 
function  of  real-time  during-  the  decElerati on  maneuver.  Therefore*  a  • 
constant  rate  of  change  of  elevator  angTe.  was:  chosen- for  this  .study. 
Rate  of  change  of  elevator  angle  was:  varied"  between -0J05.  and :-l  115; 
cleg/sec.  for  the  different  maneuvers.. 

The  effects  of  thrust  were:  studied  by  making,  two:  "test* :f Tights"  for 
each  elevator  angle  schedule..  Two  thrust  values:  were: used: for  the  e 
"flights"  which  gave  thrust  to  weight,  ratios:  of  0". 021  and.  Oil  1 '.  . 

Five  "flights"  were  made  at  -GL286T  deg/sec. .  rate. of  change :of  ' 
elevator  angle  while  varying  aircraft  weight  to:  study  the: effects :of 
aircraft  weight  on  minimum  fTying;  speed..  The:  aircraft  weight.was' 
varied  from  12,359  pounds  ta  T6\352  pounds"  iir  increments: of  -"1000  ^pounds. 

"Flights"  were  made  at  representative  rates?  of-" change  of- elevator 
angle  to  study  the  effects  of  aircraft  moment  of"  inertia  about" the 
lateral  axis,  on  minimum  flying  speed.  Values:  25".  percent  above  and : 
below  the  true  aircraft  moment  of  inertia  were:  evaluated. 

Data  from  the  original  "test  flight"  indicated  the  possible  exist- 
ance  of  an  optimum  elevator  schedule  in  determining, minimum  flying: 
speed.  An  exponential  elevator  schedule  was.  studied  and  compared  ta 
the  linear  elevator  schedule.  The  aircraft  was  decelerated  using  a 
linear  schedule  until  the  deceleration  rate  fell  below  a  preselected 
rate.  At  this  point  the  elevator  schedule  was.  changed  to  an  exponential 
schedule  until  minimum  flying  speed  was  obtained. 
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D.  DATA-ANALYSIS  METHODS 

Five  data-analysis  methods;  were:  used:  tec  determine:  mi nmum  flying, 
speed  from  the  two  main  sets  af  "fTight  tes±^'  data:..    These:methods.were: 

1 .  Constant  Deceleration  Rate 

The  linear  elevator  angle.-  schedule  resulted:  in'  near  constant : 

deceleration  rates  by  the  test  aircraft..    Minimum  flying: speed :was: 

defined  as  the  actual  minimum  speed  obtained:  during: the  deceleration. 

For  consistency,  a  deceleration"  rate  was?  defined"  as:  the; slope: of  "a: 

straight  line  drawn  from  V   »    to  T..T.  V_.    ..    This,  deceleration  rate  :was 
3  nrrn  mm 

used  in  both  the  Cr   and  C/'  methods..    Plots  were  made:  of- mini  mum  flying 
speed  as  a  function  of  deceleration  rate  and:  minimum  flying:  speed  .as  .a 
function  of  rate  of  change  of  elevator  angle  for  the:  two.  values : of  ■' 
thrust. 

The  constant  deceTerattcn  rate  method:  was:  used: to. analyze  data 
when  the  effects  of  aircraft  weight  and  moment  of"  inertia  on  minimum 
flying  speed  were  studied. 

2.  C£  Methods 


C/  was  defined  as  the  aircraft  left  coefficient: independent: 
of  normal  flight  path  acceleration: 

W 
L     qS 

The  assumption  of  aircraft  weight  and  aircraft  lift  being  equivalent 

throughout  the  maneuver  corresponds  to  considering  Cj*  as  being  an 

apparent  lift  coefficient  rather  than  a  true  value.  Minimum  flying 

speed  was  defined  by  the  C/  method  as  the  speed  or  the.  aircraft  when 

C/  was  at  a  maximum  value..  PTots  were  made  of  C.1    as:a:function  of 

max 


2T 


deceleration  rate,.  £"'   as:  a  function'  of-' deceleration  rate :and  :  C. 

max  mi  n 

as  a  function  of  rate  of  change-  ccfed'eyator  angle. 

3.  C,    Method 

(L  was  defined:  as  being;  the-  true:  aircraft  lift.coefficient  which 
compensated  for  effects  of"  normal'  flight. path  acceleration: 

W,a: 

Minimum  flying  speed  was  defined:  by  the:C,   method  as:the:speed  of  the 

r  ' 
aircraft  when  (L    was  at  a  maximum  value..    Plots. were  made  of    L"        as 
L  max 

a  function  of  deceleration'  rate^.    C,   .     as:  a  function  of  "deceleration 

mm 

rate  and    C   ..    as  a  function  of  rate  of  change  of  elevator  angle. 

mm  3  3 

4.  ATtitude  Break  Method 

The  aircraft  altitude,  trace:  had"  a  relatively  linear  character- 
istics sTope  during  the  deceleration  maneuver.  With  the:onsetof  stall 
the  characteristic  slope  would  increase  significantly  in  value.  The 
minimum  flying  speed  was  defined  as:  the  speed  of  the. aircraft  where  the 
linear  characteristic  slope  of  the:  altitude  trace  could  no  longer  be 
maintained.  The  characteristic  slope  increase  was  indicative  of 
increased  altitude  loss. 

A  deceleration  rate  was  defined  as  the  slope  of  a  straight  line 

drawn  from  h.  to  1.1  h  .  .  Plots  were  made  of  hm,-  as  a  function 
mm        mm  mm 

of  rate  of  change  of  elevator  angle.. 

5.  1-g  Break  Method 

The  minimum  flying  speed  defined  by  this  method  was  the  last 
speed  at  which  the  aircraft:  could:  maintain  1-g  flight;  This  method 
proved  unreliable  in  data  reduction  for  an  aircraft  trimmed  in  steady, 
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descending  flight  and  was  thus  not  pursued:  further..    Far- this:  method, 
an  aircraft  must  be  decelerated  from  an  fnttral"  climbr  condition.. 
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ML.     RESULTS  AND  DISCUSSION 

A.     ABSOLUTE  MINIMUM-  FLYING  SPEED! 

1.     Constant  Deceleration  Rate'  Method' 

The  minimum  flying  speeds  determined: were. plotted: in  Figure  3 
as  a  function  of  deceleration  far  thrust- to -weight:  ratios,:  0.02  .and  0.11 . 
For  a  thrust-to-welght  ratio  af  OC.OZ  the:  absolute  minimum  flying  speed 
was  83.7  knots  and  occurred  at  ct  deceleration  rate: ofc2160: knots/sec. 
A  thrust-to-welght  ratio  of  Q..TT  decreased,  the: absolute rminimura  flying 
speed  to  8G.6  knots  with  ncr  appreciable  change  in  deceleration- rate. 
Figure  2  aXso  Indicated  ab salute  minimum  fiying;  speed:was:sensitive  to 
variation  in  optimum  deceleration  rate..    The:  addition  of -"thrust  reduced 
this  sensitivity  somewhat.. 

Figure  4  indicated  absolute  minimum  flying; speed. was .relatively 
insensitive  to  rate  of  change  of  elevator..    Ratesbetween  -Q.7~and  -1.1 
deg/sec.  were  able  to  give  a  close:  approximation- of  "absolute  minimum 
flying  speed. 
•2.     C£  Method 

Maximum  values  of  C.    were  plotted  in  Figure.  5.as  a  function  of 

deceleration  rate.      C .     was  plotted  as  a  function  of "deceleration  rate 

mm    v 

in  Figure  6  and  as  a  function  of  rate  of  change  of  elevator  angle  in 
Figure  7.  A  C.'  of  1..815  was  the  maximum  attainable  in  steady, 
unaccelerated  flight.  The  absolute  minimum  flying  speed  was  83.5  knots. 
The  thrust-to-weight  ratio  of  OT.TT  decreased:  the  absolute  minimum  flying 
speed  three  knots  with  minimal"  effects,  on  deceleration  rate.   C*.  was 
fairly  sensitive  to  deceleration  rate  and.  insensitive  to. rate  of  change 
of  elevator  angTe. 
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Figure  3.     Minimum  Flying  Speed  as  a  Function  of 
Deceleration  Rate. 
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3.  CL  Method 

Absolute  minimum  flying  speed  was  88  knots:  far  a:  thrusirto:-- 

weight  ratio  of  0.02,  Figures  8  and  9..  Increasing  thrush  decreased: 

Y 
this  speed  2.5  knots.   C  .was  sensitive  tec  both  deceleration  rate:  and: 

rate  of  change  of  elevator  angle.  Figure  TO  indicated:  thee  abso:l"ute: 

maximum  C,  value  was  insensitive  to  deceleration  rata  and:  increased: 

thrust. 

4.  Altitude  Break  Method 

V 

h  i     was  plotted  as  a  function  of  deceleration-  rate:  in  Figure: 

11  and  as  a  function  of  rate  of  change  of  eJevatar  angle:  in-  Figure:  121 

The  absolute  minimum  flying  speed  for  a  tJirusrt-ta-wetght.  ratio:  of-" 01 02: 

was  83.5  knots  at  a  deceleration  rate  of -£..6  knots/sec:.     Figure:  11" 

indicated  that  two  values  for    h   •     existed  for  one  value:  of  "deceleration 

mm 

rate  between -2.6  and -3.76  kts/sec  Figure  T£  indicated:  h  •  was: 
insensitive  to  rate  of  change  of  elevator  angle.. 

The  altitude  break  method  was  dependent  on  graphical"  interpreta- 
tion for  data.  This  accounts  for  the  smaTT  amount  of" data:  scatter  in 
Figures  11  and  12. 

5.  Comparison  of  Methods 

The  general  relationships  of  absolute  minimum  flying  speed: as 
determined  by  the  various  methods  are  indicated  in  Figure  13.. 

The  constant  deceleration,  C.  and  altitude  break  methods  defined 
the  absolute  minimum  flying  speed  to  be  83  knots,  Figure  13.  The  C. 
method  absolute  minimum  flying  speed  was  88  knots.  Increased  thrust 
decreased  the  absolute  minimum  flying  speed  defined  by  the  four  methods. 
about  three  knots. 
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Figure  8.     Minimum  Flying  Speed  as  a  Function  of  Deceleration 
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Figure  13.  Comparison  of  Aircraft  Parameters  at  Absolute 
Minimum  Flying  Speed  T/W  <  0.02. 
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The  deceleration  rate  was  approximately  -2L6i  knots/see:at:the= 
absolute  minimum  flying  speed  far  all  four  methods:. .  Sjgpificant:changes 
in  deceleration  rate  due  to  increased  thrust  was;  observed: only  irrthe. 
altitude  break  method. 

Minimum  flying  speed  determined  by  the  C.  method: occurred  :at: 
the  smallest  rate  of  change  of  elevator  angle*.  -O'-^  degs/see: .  Minimum 
flying  speed  determined  by  the  constant,  deceleration,,  C.1  and:altitude: 
break  methods  occurred  at  -Q..9  deg/ sec.  rate  ccf"  change, of" elevator  angle. 
Absolute  minimum  flying  speed  defined  by  the  C  method. occurred. at :the: 
smallest  elevator  angle,  approximately  -9  degrees:,  whereas:  approximately 
-13  degrees  elevator  angTe  was  required  for  the:  other  methods: . 

The  angle  of  attack  at  the  absolute  minimum  flying;  speed  .was 
approximately  28  degrees  for  the  constant  deceleration, .  Cj*  and.altitude 
break  methods,  and  considerably  Tess  at  ZT  degrees  for  the: C.  method. . 
Thrust  effects  on  angle  of  attack  at  the  absolute:  minimum  flying. speed: 
were  insignificant  (less  than  Q.5  degrees)  for  all  methods. 

Elevator  angles,  lift  coeffT'cients:  and  angle:  of' attacks: at: 
minimum  flying  speed/stall  speed  were  less  than,,  in  some  cases;  the: 
maximum  static  values  defined  for  the  aircraft.  Possible  explanations  . 
are: 

1.  Rate  of  change  of  elevator  angle  was  programmed: to  become 

zero  two  seconds  after  V  •  was  obtained.  For  small  values  of  elevator 

rate,  V  •  occurred  much  sooner  than  the  maximum  obtained  values: of  C,  , 

o  and  6  .  Therefore,  elevator  angle  was  held  constant  until  the  end  of 

the  maneuver,  restricting  the  aircraft  from  obtaining  the  maximum  values 

of  Ci  »  a  and  5  . 
L        e 
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2.  The  variations  of  drag  coeff f event:  with  angl'ee of-  attack-  used: 
in  the  table  look  up  were  representative  values:  and:  not: necessarily  true 
values  for  the  F-94A.  Reduced  drag  coefficients:  would:  there  fore  effect: 
the  occurrence  of  minimum  flying  speed. 

3.  In  defense  of  the  Tow  values  of  C  ,.  a  and:  Si..  it:must:be: 

u.  e. 

remembered  that  this  is  a  dynamic  study  of  stalT  traits:,.  Actual"  flight: 
test  data,  [Ref.  4],  indicates  angle  of  attack  is:  less:  than  static: stall 
angle  of  attack. 

B.  FAR  MINIMUM  FLYING  SPEED 

The  minimum  flying  speed  defined  by  Federal  Aviation  Regulations: 
corresponds  to  the  minimum  speed  obtained  at  a  one  knot/sec:.,  decelera- 
tion rate.  The  FAR  minimum  flying  speed  determined  by  the  constant: 
deceleration,  C.  and  C  methods  was  approximately  105  knots.  The: 
altitude  break  method's  FAR  minimum  flying  speed,  was:  1D9.:  knots. .  Increased 
thrust  decreased  the  altitude  break  minimum  speed  two  knots  and  decreased 
the  minimum  speed  defined  by  the  remaining  methods  less. than  one: knot. 

The  rate  of  change  of  elevator  was  -  O..0S  dag/ sec-,  at- minimum  flying 
speed  determined  by  the  constant  deceleration  rate,.  C.  and_ C^  methods. 
The  altitude  break  method  rate  of  change  of  elevator  angle  was  -0.05 
deg/sec. 

Angle  of  attack  was  approximately  15  degrees  and  elevator  angle  -4.8 
degrees  for  all  methods. 

C.  COMPARISON  OF  MINIMUM  FLYING  SPEEDS 

The  absolute  minimum  flying  speeds  were  categorically  20  knots  less 
than  those  defined  by  FAR  one  knot/sec.  deceleration  rate,  Figure  14. 
Thrust  had  more  effect  on  the  absolute  minimum  flying  speed  whereas  rate 
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Figure  14.     Comparison  of  Aircraft  Parameters  at  FAR  Minimum 
Flying  Speed. 
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of  change  of  elevator  angle  had  mare  effect  an  the  FAR  minimum  fTying: 
speeds.  Angle  of  attack  and  elevator  angTe  were  significantly  greater 
with  the  absolute  minimum  flying  speeds. 

D.  SENSITIVITY  STUDIES 
1 .  Aircraft  Weight 

For  a  representative  rate  of  change  of  eTevator  ang^Te,.  Figure: 
15  indicates  minimum  flying  speed  was  a  linear  function  of"  aircraft-,, 
weight. 

From  the  definition: 

W  . 
C, 


hnax.     q  ■  S 

minimum  flying  speed  is  seen  to  be  a  function  of  the  one-haTf  power- of- 
weight.  Actual  flight  test  data  indicates  minimum  flying  speed:  is:  a: 
linear  function  of  aircraft  weight  [Ref.  4]. 

2.  Aircraft  Pitch  Moment  of  Inertia 

When  the  aircraft  pitch  moment  of  inertia  was:  increased  or 
decreased  25  percent,  there  were  no  significant  effects-  on  minimum 
flying  speeds  or  other  aircraft  parameters.  This  conclusion  was  not- 
based  on  a  systematic  study  but  rather  a  brief  look  at  three  data 
runs.  The  trends  were  present  but  a  complete  study  is  warranted.. 

3.  Exponential  Elevator  Rate 

Comparative  plots  of  linear  and  Tinear  plus  exponential  elevator 
rates  were  made  in  Figure  16  for  initial  Tinear  rates  of  -0.06  and 
-0.115  deg/sec.  The  exponential  rate  was  ex  where  x  was  increased  0.04 
each  second.  V  .  obtained  by  the  linear  pTus  exponential"  eTevator  rate: 
was  significantly  less  than  that  obtained  by  the  Tinear  elevator  rate. 
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These  results:,,  thcrugir    arrVy.  et  sampling  of'' data : , .  tend  : to  support  the, 
existence  of  an  optimum"  eTeyaetor'  schedule  in  determining  :minimum  flying 
speed. 
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IV.     CONCLUSIONS  ATlDt  RECDMMEND'ATIDNS: 

Digital  simulation  has  great,  potential"  in- the:  evaluation  of 'minimum 
flying  speed.     The  options  available  are  practicably  unlimited. 

The  determination  of  absolute  minimum  flying:  speed: and: aircraft: 
behavior  with  digital  sfmuTatfarr  will"  enable  the:  engineer  to:  better 
predict  the  aircraft  charactertsttcs  encountered  during; test:f light; 
The  test  pilot's  job  will  be  snimpTified  somewhat,  if" he: has: knowledge  of 
general  aircraft  characteristics  prior"  to  flight.. 

The  operational  pilot's  performance  should  be:  enhanced: by  having 
an  insight  into  his  aircraft's  performance  beyond:  the:  certified  stall 
speed. 

The  computer  program  coded  for  this  study  may  be: used. for  many. 
fields  of  study  as  indicated  by  Ref..  7\:  where  it:  was:  used:  to:make:a: 
stability  analysis  in  the  vicinity  of  minimum  flying  speed. 
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AFFEflCUX  K 
NON-LINEAR  EQUATIONS'  OF  AIRCRAFT"  MOTION' 

I.  Euler's  Equations  of  Aircraft  Motion. 

The  aircraft  was  considered  as:  &  rigid: body  with  a- Cartesian 
coordinate  system  fixed  at  the  center-  of-'  mass:,-  Figures  17"..  The -frame  of 
reference  Cxyz,  a  body  a*is  system,,  moved  with  the:  aircraft:and:Cx  was 
fixed  to  the  longitudinal;  axis  erf  the  aircraft-.    TTie: resultant  .external 
forces  and  the  moments  of  these:  forces:  about:  the  aircraft:center  of  mass, 
when  referred  to  the  frame  of  reference:  Cxyz:,.  where:    tbe:vector  equations 
of  motion: 

m         _    _ 

F  =  m  tt—  +  m  ox  x  V 

sr=§  ^»r 

The  scalar  components  of  these  equations  are:  Euler's: equations: of  ~ 
aircraft  motion. 

II.  Aircraft  Orientation  and  Axes. 

Reference  frame  Cxyz  could  not  be  used  to  describe. the  position  and 
orientation  of  the  aircraft  as  it  was  fixed  to  the  center  of  mass  and 
moved  with  the  aircraft.  An  earth  fixed  reference  frame,  assumming 
the  earth's  rotation  was  negligible,  was:  chosen  to  describe  the  aircraft 
motion. 

Euler's  equations  of  motion  are  valid  for  any  orthogonal  reference 
frame  fixed  to  the  aircraft  center  of  mass:..  To:  simplify  the  .equations 
of  motion  and  expressions  for  aerodynamic:  forces,. the  reference  frame, 
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Cx'y'z',  was  chosen,  as  stab  fifty,  axes:-,.  f  ..e£..,.  Go'   pointed:  iir  tbeedireetion- 
of  motion  of  the  aircraft  center  ccf  mass,.  Ff-gure:  IT.. 

III.  Assumptions 

In  addition  to  the  assumptions  previously  made:- defining: aircraft: 
orientation  and  axes  ft  was  assumed,  that.  a:ircra.ft- moti on  was:  restricted: 
to  the  Cxy,  longitudinal  plane..  Thfs  assiimptiirrr  equated:  the  following: 
quantities  to  zero:  L,  N,  F,.  Rr  $:,.  «„  V.  and:  Y..   Cxy  was;  also: assumed  to: 
be  a  plane  of  symmetry  which  equated  the  products of"  inertia^. Dland:F^ . 
to  zero. 

IV.  Non-Dimensionalfzation 

The  external  aerodynamic  forces?  and.  moments:  were:  expressed -as: 
functions  of  dimensionless  force  and  moment  ctxeffici en ts.,.  dynamic; pressure, 
aircraft  planform  wing  area  and  dt  reference:  length..  These: equations, . 
in  scalar  form,  were: 

When  these  values  were  substituted  into  the:  scalar  resultant  external 
force  and  moment  equations  the  following  equations  were  obtained: 

p  Vc  S  C* 

II  =  = —  -  5  sin  &  -  QW 

2m 

P  I?  S  C_ 

W  = =-  +  g  cos  e  +  QU 

2m 
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BODY  AXES 
STABILITY  AXES 


Figure  T7.  Aircraft  Axes 
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p  V?  S  c  C 

q  =  S S- 

B 

q  =  a 

V.  Aerodynamic  Force  and  Moment  Coefficients: 
From  Figure  18,  the  force  coeff ici  ents:  were: 

C       =    Cj-    +     (I    ssi'tv  x  —  C~  cos:  a: 

Cz     -  -C   cos  a  -  C.  srin     cc 

The  pitching  moment  and  Ttft  coefficients:  consisted- of-" linear  portions. 
The  non-linearity  was  present  in  C,(a)  and  Cm(a) ,. which  were.non-linear 
curves  of  C,    and  C     versus     ..    Complete  expressions,  for  the: total " 
pitching  moment  and  lift  coefficients  were: 


f[ 

1  Cm    e:  -<-  r    . .  a. 
m              rrr 

i      Q                 a:.    J 

1 

m  -  .       e  t 
J  •      6ee 

L    =    V^     +CL^    *  +  (lw6*.    ' 


VI.  Final  Equations 

Angle  of  attack  appeared  in  the  pitching  moment  equation  as  a 
derivative.  To  obtain  a  set  of  simultaneous  differential  equations  that 
could  be  solved  by  one  of  the  methods  available,  another  equation  was 
added  defining  the  rate  of  change  of  angle  of  attack.  The  three-degrees- 
of-freedom,  non-linear  differential  equations  of  aircraft  motion  programmed 
for  the  Naval  Postgraduate  School"  IBM  360  digital  computer  were: 
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Figure  18.     Aerodynamic  Force  and  Moment  Coefficients 
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u    = 


p  *?  s 

c 


2m 

g  sin  6-  -  QW 


C^-HU^a)    -H  (T     en-CL        6U  isin  a- -Cp 


COS   a 


6£ 


W     = 


^VcS    I 


Znr      | 
+  §  cos  8F ■  +  QU 


fc(a) 


C.       6y  jcos:ai+-CD  sin  a 


Jl 


=    Q 


Q    = 


2B 


C   (a)   +.-  oTr |C  .    a 
V      a 


m.  .     e 

6ee 


W  cos  a  -  U  sin  a 
a     =     
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APPENDIX  BI 
DEVELOPMENT  0Fr  COMPUTER:  PROGRAM 

I.     COMPUTER  FACILITY 

The  original"  intent,  was,  to:  program  tbe:equations  :of  motion  on  the 
Aeronautical  Engineering;  department. EAL 580! analog  computer.     Problems 
with;  insufficient  equipment"  toe  program  theiequations  :of  motion  on  the 
analog  computer  necessitated  the:  use:of-digitarsimulation. 

The  Naval  Postgraduate  School's: IBM  360! digital  computer  offered 
mar.y  advantages  and:  was:  chosen  as:  the:  research  tool;     Several  major 
advantages  were::  the.  Targe:  amount; of"  core-:avai Table,  .dimensioned 
quantities  could  he  used:,,  plot  subroutines^were  available  for  time 
Biistory  studies  and  numerous  methods  werecavailable  .for  solving 
stimuli  taneaus  non-Tinear  differential'  equations  on  a  :digital  computer. 

IE.  COMPUTER  PROGRAM 

Of  the  many  one:  step  procedures:  that:  were  -.available  :for  solving 
systems  of  simultaneous  firs.t-order  ordinary  differential  equations, 
a  fourth-order  Runge-Kutta  algorithm  was: chosen,  [Ref.'  8}."  This  method 
offered  simplicity,,  accuracy  and:  short" computer  time.  The  Runge-Kutta 
procedures  were  programmed  into  subroutine  RUN  KITH  RUN KUT  called 
function  subroutines  Fl ,  F2,  F3,  F4  and  F5  to  calculate  the  derivative 
terms  at  time  T(I),  T(I")-t-HI/2  and  T(I)+HI.  The. derivative  terms  were 
multiplied  by  the  appropriate,  weight: factors  and  added  to  the  variables 
calculated  at  time  T(I)  which  gave  values,  for  the  variables  at  time 

t(;ih);- 
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Function  subroutine  FT,.  caTTet  saitrroiitine:  SPUN  to:calculate:the 
non-1  inearities  fn  C.*  C    arret  Cy..    SPUN  was:  obtained:ffom' the  scientific 
subroutine  package  of  the  NavaT  Postgraduate:- School'  IM  36Q:digitat" 
computer  and  modfffed  for  this  program.     Non-Linear  tables :of~'Cp,  .Cffl  and 
CD  as  a  function  of  angTe  of  attack  were  read:  into:  SPUN.  .   When  supplied 
values  of  angle  of  attack,.  SPLIN  used  a:  cubic: spline: function  to  provide 
interpolated  values  of  (L  ,.  C    and  £>.. 

Initial  aircraft  trim  conditions'  U,.  14,.  e,.  §"».  a,  and:6.  were:supplied 
to  RUNKUT  at  time  zero  for  initiation  of  the:  step:  procedure.  .  These 
values  were  calculated  in  Subroutine:  TRIM  from  stati c: force :and  :moraent 
equations. 

Control  of  the  aircraft  was  accomplished: through  changes :in  elevator 
angle  and  thrust  which  simulated  control"  stick  and  throttle. movements. 
The  options  available     were  seTecterL  through  the: integers:on- 
the  first  data  card. 

III.     TRUNCATION  ERROR  AND  STEP-SIZE  ANALYSIS. 

Error  analysis  by  rigorous  mathematical  derivationssaFe:virtualIy 
impossible  to  implement  for  higher-order  Runge-Kutta  algorithms  .for 
systems  of  differential  equations.     The  truncation  error: 

-  II    l  \ 

et  "  15     lyn+l,2  ~"yn+l,l; 

for  a  first-order  ordinary  differential  equation  was  applied  to  the 

system  of  differential  equations  as  a  step-size  control  mechanism,  [Ref.  8] 

Step-size  and  truncation  errors  are  listed  in  Table  II. 
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Table  EI 


Local  Truncation  Errors 


Step-Size  Error 

0.2  -2?. 521  x  m"^ 

0.1  -3EIS57X  ltf€: 

0.05  11.56'  *  iu~  ^: 

0.025  -33.SE    x>  1D:"& 
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13.  ABSTRACT 


Aircraft  minimum  flying  speed,  as  determined  by  actual  flight  test, 
is  published  in  aircraft  handbooks  for  pilot  guidance.  The  test  flight, 
results  are  used  to  determine  and  confirm  take-off  and  landing  speeds., 
field  lengths,  left-hand  portion  of  the  maneuvering  envelopes  (V-n  diagram), 
etc.  Determination  of  the  absolute  minimum  flying  speed  of  an  aircraft 
on  the  other  hand,  has  not  been  of  prime  importance  in  flight  test. 

In  the  present  analysis  digital  simulation  allowed  the  systematic 
study  of  not  only  the  minimum  flying  speed  as  defined  by  Federal  Aviation 
Regulations  but  also  the  absolute  minimum  flying  speed  attainable  in  steady, 
unaccelerated  flight.  The  study  included  such  effects  as  deceleration 
rate,  rate  of  change  of  elevator  angle,  aircraft  weight  and  pitch  moment 
of  inertia. 

It  was  found  for  an  assumed  light-weight  fighter  aircraft  that  the 
absolute  minimum  flying  speed  was  approximately  2Q  knots  less  than  the  FAR 
minimum  flying  speed.  Moreover  the  FAR  minimum  flying  speeds  tended  to  be 
quite  sensitive  to  rate  of  change  of  elevator  angle. 
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